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Abstract 

Quaternion-valued wireless communication systems have been studied in the past. Although progress has been made 
in this promising area, a crucial missing link is lack of effective and efficient quaternion-valued signal processing 
algorithms for channel equalization and beamforming. With most recent developments in quaternion-valued signal 
processing, in this work, we fill the gap to solve the problem by studying two quaternion-valued adaptive algorithms: 
one is the reference signal based quaternion-valued least mean square (QLMS) algorithm and the other one is the 
quaternion-valued constant modulus algorithm (QCMA). The quaternion-valued Wiener solution for possible block- 
based calculation is also derived. Simulation results are provided to show the working of the system. 

Keywords: Polarisation diversity, four-dimensional modulation, quaternion valued signal processing, channel 
equalization, beamforming, constant modulus, least mean square. 


1. Introduction 

Increasing the capacity of a wireless communication system has always been a focus of the wireless communi¬ 
cations research community. It is well-known that polarisation diversity can be exploited to mitigate the multipath 
effect to maintain a reliable communication link with an acceptable quality of service (QoS), where a pair of anten¬ 
nas with orthogonal polaristion directions is employed at both the transmitter and the receiver sides. However, the 
traditional diversity scheme aims to achieve a single reliable channel link between the transmitter and the receiver, 
while the same information is transmitted at the same frequency but with different polarisations, i.e. two channels. 
This is not an effective use of the precious spectrum resources as the two channels could be used to transmit different 
data streams simultaneously. For example, we can design a four-dimensional (4-D) modulation scheme across the 
two polarisation diversity channels using a quaternion-valued representation, as proposed in ll]]. An earlier version 
of quaternion-valued 4-D modulation scheme based on two different frequencies was proposed in However, due 
to the change of polarisation of the transmitted radio frequency signals during the complicated propagation process 
including multipath, reflection, refraction, etc, interference will be caused to each other at the two differently polarised 
receiving antennas. To solve the problem, efficient signal processing methods and algorithms for channel equalization 
and interference suppression/beamforming are needed for practical implementation of the proposed 4-D modulation 
scheme. 

Recently, quaternion-valued signal processing has been introduced and studied in details to solve problems related 
to three or four-dimensional signals such as vector-sensor array signal processing flElSHSIH], and wind profile 
prediction ilii. With most recent developments in this area, especially the derivation of quaternion-valued gradient 
operators and the quaternion-valued least mean square (QLMS) algorithm |[i3,[ni[il, we are now ready to effectively 
solve the 4-D equalisation and interference suppression/beamforming problem associated with the proposed 4-D mod¬ 
ulation scheme. Now the dual-channel effect on the transmitted signal can be modeled by a quaternion-valued infinite 
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Figure 1: Wireless communication employing polarisation diversity, where both the transmitter and the receiver sides are equipped with a pair of 
antennas with different polarisation directions. 


impulse response (HR) or finite impulse response (FIR) filter. At the receiver side, for channel equalisation, we can 
employ a quaternion-valued adaptive algorithm to recover the original 4-D signal, which inherently also performs an 
interference suppression operation to separate the original two 2-D signals. Moreover, multiple antenna pairs can be 
employed at the receiver side to perform the traditional beamforming task to suppress other interfering signals. 

In particular, two representative quaternion-valued equalisation/beamforming algorithms will be derived; the first 
one is the quaternion-valued Wiener filter as a follow-up to the previously derived QLMS algorithm for reference 
signal based equalisation/beamforming, and the second one is the quaternion-valued constant modulus algorithm 
(QCMA) for blind equalisation/beamforming. Compared to the summary contribution in iflUl . in addition to the 
detailed analytical modeling steps, the main difference is the GCMA algorithm and the related simulations. Although 
quaternion-valued wireless communication employing multiple antennas has been studied before, such as the design 
of orthogonal space-time-polarization block code in Cl, to our best knowledge, it is the first time to study the 
quaternion-valued equalization and interference suppression/beamforming problem in this context. Moreover, the 
dual-polarised antenna pair or an array of them has a similar structure to the well-studied vector sensors or sensor 
arrays |T1 H 12. iSl, where they are used mainly for traditional array signal processing applications. Although the 
recently developed quaternion-valued array simal mocessing algorithms based on such traditional array applications 
employed a quaternion-valued array model l^lsllalTI], the desired signals are still traditional complex-valued signals, 
instead of quaternion-valued communication signals. 

In the following, the 4-D modulation scheme based on two orthogonally polarised antennas will be introduced in 
Sec. |2]and the required quaternion-valued equalisation and inter-channel interference suppression solution and their 
extension to multiple dual-polarised antennas are presented in Sec. |3] Simulation results are provided in Sec. 01 
followed by conclusions in Sec. |5] 


2. Quaternion-Valued 4-D Modulation 


In traditional polarisation diversity scheme, as shown in Fig. [T] each side is equipped with two antennas with 
orthogonal polarisation directions and the signal being transmitted is two-dimensional, i.e. complex-valued with one 
real part and one imaginary part. In the quaternion-valued modulation scheme, the signal is modulated across the 
two antennas to generate a 4-D modulated signal. Such a signal can be conveniently represented mathematically by a 
quaternion il9i|2C 


. A quaternion is a hypercomplex number defined as 
q^qQ + iqi + jqi + kqj, , 


( 1 ) 


where qo is the real part of the quaternion, and qi, q 2 and q^ are the three imaginary components with their corre¬ 
sponding imaginary units /, j and k, respectively. The conjugate of a quaternion, denoted by q*, is defined as 

q* ^qo- iqi - jq 2 - kqj, . (2) 

/, j and k satisfy the following conditions 


a ^ jj ^ kk ^ -I , (3) 

ij = -ji ,jk^ -kj = i , ki = -ik = j . (4) 

As a result, quaternionic multiplications are noncommutative. 

As an example, corresponding to the 4-QAM (Quadrature Amplitude Modulation) in the two-dimensional case, 
for the 4-D modulation scheme, qo, qi, q 2 and qT, can take values of either 1 or -1, representing 16 different symbols. 
We can call this scheme 16-QQAM (Quaternion-valued QAM) or 16-Q^AM. 
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Figure 2: Channel model and the reference signal based equalizer for quaternion-valued signals. 


3. Quaternion-Valued Equalization and Interference Suppression/Beamforming 

3.1. Channel model 

The signal transmitted by the two antennas will go through the channel with all kinds of effects and arrive at 
the receiver side, where the two antennas with orthogonal polarisation directions (Note that orthogonal polarisation 
may not give the best performance for a specific scenario) will pick up the two signals. Again the four components 
of the received signal can be represented by another quaternion. We use s,[n] and Sr[n] to represent the transmitted 
and received 4-D quaternion-valued signals, respectively. Then the channel effect can be modeled by a filter with 
quaternion-valued impulse response fc[n], i.e. 

Sr[n] = s,[n] * fc[n] + qa[n] , (5) 

where qa[n] is the quaternion-valued additive noise, as shown in Fig.|2] 


3.2. Reference signal based approach 

To recover sfn] from Sr[n] or estimate the channel, as in the 2-D case (complex-valued), we can design a 
quaternion-valued equalizer. One choice is a reference signal based equalizer, among many others corresponding 
to the complex-valued case. Now assume we have a reference signal r[n] available. Then we can employ the standard 
adaptive filtering structure shown in the second half of Fig. and update the equalizer coefficient vector w with a 
length of L by minimising the mean square value of the error signal e[n] imoEi. 

The cost function is given by 

J[n] = e[n]e*[n] , (6) 

where 

e[n] - r[n] -y[n] - r[n] - w^Sr[n] . (7) 

with w being the equalizer coefficient vector and Sr[n] holding the corresponding received signal samples from Sr[n] 


S/'[n] = [Sr[n], Sr[n - , Sr[n - L + \]]^ 


( 8 ) 


For a general quaternion-valued function f{q) of a quaternion q, the gradient of / with respect to q is defined as 


df 

1/5/ 

df. 

df . 

df A 

dq 

4 

dqi 

dqi^ 

dq3 ) 


dq* 


^ + Ki. 

dqo dqi 


dqi 


(9) 

( 10 ) 


Following the derivations in we have the gradient of J{n] with respect to the coefficient vector as 

follows 

Vw7[n] =-is,[n]e*[n], (11) 
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which leads to the following update equation for the coefficient vector with a step size of yu, i.e. the QLMS algorithm: 


w[n+l] = w[n]-MVw^W)* (12) 

= y/[n\+^i{e[nX[n]). (13) 

For a solution equivalent to the classic Wiener filter in the complex-valued case, based on the instantaneous 
gradient result of (HB, the optimum solution should satisfy 

E{-\sAn]e*[n]} = 0, (14) 

i.e. 

E{sAnV[n]] = £{s,[n]r*W-s,[n]sf[n]w:^,) 

= p-R.Xpr^O, 

where the cross-correlation vector p = £’{Sr[n]r*[n]} and the covariance matrix X = £{Sr[«]sf [«]} 

^swipt = p => <pi = r;/p■ 

We can use the above equation to obtain the optimum weight vector directly. 

Now consider the multiplication of two quaternions a and b with c - ab. They are expressed as 

a — uq + ia\ + ja 2 + ka^, 

b — bo + ib\ + jb2 + kb^ 

c - Co + ici + jc 2 + kco (17) 

where the subscripts indicate the corresponding components of the quaternion. Then we have 


(15) 

. Then we have 

(16) 


Co 

= booo - biai - b2a2 - b^aj. 


Cl 

- booi + b\aQ - b2a2 + b2a2 


C2 

- boa2 + biao, + b2ao - b^ai 


C3 

- boa^ - b\a2 + b2a\ + b^ao 

(18) 


Using this result, we can obtain the solution to (fTST l using real-valued matrix operations. First we define the following 
two vectors 


w, 


opt 


P = 


WKpuo)^ -JiKpoo)^ -K{Kp,,o)^ 

■ • ■ ^Ri<pUL-y)^ -JiKpOL-^ -^«p,L-l)]l^ 

[/^(p)^/(p)^7(p)^/:(p)^]^ 


(19) 


where Wopi,i, 
operation of 
also define 


with 


/ = 0, • • ■ , L - 1 is the /-th element of the optimum weight vector Wop?, R{-), /(•), -/(•), and K{-) are the 
taking the real and three imaginary components of the quaternion inside the brackets, respectively. We 


Rt../ — 


Ri, - [R.s„o; • ■ • 



(20) 

Riri) -/(r,) 

-Jiri) 

-^(r,) ) 


/(ri) 

R{ri) 

-Kin) 

+Jiri) 

(21) 

m) 

Kiri) 

Riri) 

-m 

Kiri) -Jiri) 

Kn) 

R{ri) , 



for 1 = 0, l,-'-,L-l, where r/ is the 1-th column vector of the covariance matrix R,,.. 
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Then, according to (fTSl l. we have the following relationship 


Rj.Wopr = P , 

where all the matrix and vectors involved are real-valued. Then Wopt is obtained by 

Wop, = R,/p . 

From Wopt, we can then easily deduce Wop,- 


(22) 


(23) 


3.3. Constant modulus based approach 

When a reference signal is not available, it is still possible to perform equalisation and interference suppres¬ 
sion/beamforming by employing other properties of the signals. An algorithm designed to work without knowledge 
of the transmitted signals falls into the category of blind equalisation and beamforming approaches iHllllIl. One 
representative blind equalisation algorithm in traditional communication systems is the constant modulus algorithm 
(CMA) 

There are many variations to this algorithm and the basic form is based on minimizing the following cost function 




(24) 


where y[n] is the output of the equalizer and y is the dispersion constant, defined by y = being symbols 

of the modulation scheme. 

For our quaternion-valued 4-D wireless communication system, we can develop a similar quaternion-valued CMA 
(QCMA) for blind equalization and beamforming. Taking the gradient of Jcm with respect to the quaternion-valued 
coefficient vector w, and using the third chain rule of the restricted HR gradient operation given in 11211 as the inter¬ 
mediate function is real-valued, we have 

1 „ <5(y[n] • y*[n] — y) 

VwTcm = -£{(yW-y*W-y) ^ ^ 

1 r-.. r 1 *r 1 ^d(y[n] ■ y*[n]) 

= ^E{(y[n] -y [n] - y)---} 

2 aw 


^ £■((}'[«]•/[«] -r) 


(9(w^Sf[n]sf [n]w*)j 

5^^ * 


Using the result of ifioll . we have 


(9(w^Sr[n]sf [n]w*) 1 

dw 2 


S.[n]y*[n] 


Then we have 


VwTcm = -^E{{y{n\ ■ f[n] - y)s,[n]y*[«]} 


(25) 


(26) 


(27) 


Using the instantaneous gradient V^Tcm to replace V^Jcmj we then obtain the final update equation for our quaternion¬ 
valued constant modulus algorithm 


(28) 
(29) 

where the constant j has been absorbed into the step size p. 


w[nH-l] = w[n]-piV^JcmY 

= w[n] - p(\y[n]\^ - r)y[n]s*[n] , 
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Figure 3: Multiple antenna pairs at the receiver side, where each pair is composed of two differently polarised antennas. 
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Figure 4: A MIMO system with multiple antenna pairs at both the transmitter and the receiver sides, where each pair is composed of two differently 
polaiised antennas. 


3.4. Extension to multiple antennas and MIMO systems 

We can extend this design to multiple antenna pairs for beamforming to suppress other quaternion-valued interfer¬ 
ing signals, as shown in Fig.[3j or Fig. |4]for a general multiple-input-multiple-output (MIMO) system. The channel 
model for an M x N system is shown in Fig.|5] where s,^„,[n], m - 0,1, - ■■ ,M - 1 is the transmitted signal, while 
Srj[n], I - 0,1, • • ■ ,N- 1 is the received signal. qaj[n], / = 0,1, ■ • ■ , V - 1 is the added channel noise and fm,i[n] is 
the quaternion-valued channel impulse response between the m-th transmit antenna pair and the /-th receive antenna 
pair. 

Mathematically, we have 
where 


A 

s,[n] 

Sr[«] 

qj«] 


S;.[n] = ■ s,[n]+qa[n] , 


(30) 


/o,oW 

fi,o[n] 


fo,i[n] 

/i.iW 


/m-1,o[«] /m-1,i[«] 

[qa,o[n],qa,l[n],--- 


/m-i.n-i [«] 


(31) 


Applying an L x 1 weight vector w„,,/ = [wmj.o, Wm.ui, ■ ■ • ,Wm,i,L-iY, m = 0,1, • ■ • ,M - 1, / = 0,1, • • ■ , A - 1, to 
the received signal Sy.iin} and then combining the N corresponding outputs together, we obtain one of the estimated 
signals y„,[n], i.e. 

ym[n] = yvJjM , (32) 
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Figure 5: Channel model for a quaternion-valued MIMO system. 



Figure 6: Learning curve of the quaternion-valued equalizer using the QLMS algorithm. 


where 


Sr[«] = 

Srj[n] = [Srj[n],Srj[n- I],--- ,Srj[n- L+l]f . (33) 

The optimum weight vector for Wm can be obtained using either the QLMS algorithm or the QCMA introduced 
before so thaty,„[n] becomes a good estimate of one of the transmitted signals. 

4. Simulation Results 

In the following, we give three sets of simulation results. The first two are for the QLMS algorithm and the third 
one for the QCMA. 

The first set of simulations is based on the structure in Fig.|2]and it is a reference-signal based channel equalization 
problem. The signal transmitted is 16-Q^AM modulated and the SNR at the receiver side is 20 dB with quaternion¬ 
valued Gaussian noise. The channel impulse response fc[n] is a 4-tap quaternion-valued FIR filter with a Gaussian- 
distributed coefficients value, which is generated randomly for each run. The equalizer filter has a length of 15. The 
learning curve based on averaging 200 simulation runs using the QLMS algorithm is shown in Fig. |6] with a step size 
fj. - 0.000012, with about -10 dB error at the steady state, indicating a reasonable channel estimation result. Note 
that the optimum ju for each channel realization is different and a fixed p for all randomly generated channel impulse 
responses will lead to a much less favorable result, which is why the steady state error is relatively large. We have 
also calculated the corresponding bit error rate (BER) based on this non-optimum /j. and it is about 0.29%. 
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Figure 7: Scatter plots before and after equalization using the QLMS algorithm. 


However, if we focus on one specific channel response, such as the one given below 


-1.1696 

-0.1135 

0.0958 

-0.9583 


0.4989 

-0.3639 

-0.2515 

-0.6047 


1.0571 

-0.3791 

1.4957 

0.5183 


0.4911 

0.8653 

0.1552 

-0.8966 


(34) 


where each column gives the four components of the quaternion-valued coefficient of /c[n], and use the same value 
of yu = 0.000012, we will be able to achieve a zero BER result, as shown by the scatter plots of constellation before 
and after equalization in Fig. |7] Since we can not show the 4-D scatter plot directly, we have split it into two 2-D 
plots in Fig. |7] where the upper row is for the scatter plots before equalization, and the lower row for the plots after 
equalization. We can clearly see that the equalization operation has been successful. 

In the second set of simulations, we consider a 2 x 2 MIMO array and the two transmitted quaternion-valued 
signals have the same normalized power, with an SNR of 20 dB at the receiver side. All the other parameters are the 
same as the first one, except that now the step size has been changed to fi - 0.0000006. Fig.[8]shows the result, again 
a reasonable performance with such a fixed non-optimum step size. 

Moreover, we have used the MATFAB function ‘tic’ and ‘toe’ to calculate the running time for the two scenarios. 
Based on a Windows 7 computer with Intel Core i5-2467M CPU (1.60GHz) and 4GB memory, it took about 4.9 
seconds for each run of the first scenario and 10.6 seconds for each run of the second scenario. This time is dependent 
on the stepsize adopted for the algorithm. A larger stepsize will significantly reduce the time required for convergence. 
The reason for our current step size value (which is quite small) is to make sure the algorithm will converge for all 
channel realizations in the simulation. In practice, this step size could be normalized according to the received signal 
powers to give a much faster convergence rate, in a similar way to the normalized FMS algorithm in literature for 
complex-valued signals SB. 

Now we consider an example for the QCMA. A well-known problem with constant modulus based algorithm 
is that it may converge to local minima depending on its initial values. We find that it is very difficult to initialize 
the algorithm properly to perform both the equalization and beamforming tasks simultaneously based on the general 
MIMO structure. So here we only show a case for blind equalization using the QCMA. The setting is the same as the 
first set of simulations. Since we are using the 16-Q^AM modulation scheme, we have chosen y - 4. The stepsize is 

yU = 0.0012. 
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Figure 8: Learning curve of the quaternion-valued equalizer/MIMO beamformer using the QLMS algorithm. 



Figure 9: Learning curve of the quaternion-valued constant modulus algorithm. 


The learn curve in terms of the instantaneous cost function |(|y[n]p - 4)^ averaged over 200 runs is shown in 
Fig. |9] where we can observe that the approximately constant modulus status has been reached. The scatter plots 
of constellation before and after equalization for the specific channel realization in (l34l i are shown in Fig. [TO] which 
clearly demonstrate that the equalization operation has been successful. One note is the phase shift in the plot, which 
is normal as the blind equalization algorithm we employed here is ambiguous to an arbitrary phase shift. 

5. Conclusions 

A 4-D modulation scheme using quaternion-valued representation based on two antennas with different polarisa¬ 
tion directions has been studied for wireless communications. Although quaternion-valued wireless communication 
systems have been investigated in the past and research progress has been made in this promising area, a crucial miss¬ 
ing link is lack of effective and efficient quaternion-valued signal processing algorithms for channel equalization and 
interference suppression/beamforming. To fill this gap, two representative quaternion-valued signal processing algo¬ 
rithms have been introduced; one is the reference-signal based quaternion-valued least mean square algorithm with 
the associated Wiener solution, and other one is the blind quaternion-valued constant modulus algorithm. Three sets 
of simulation results were provided showing that both algorithms work effectively in the single-input-single-output 
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Figure 10: Scatter plots before and after equalization using the quaternion-valued constant modulus algorithm. 


and multiple-input-multiple-output cases and such a 4-D modulation scheme can therefore be considered as a viable 
approach for future wireless communication systems. 
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